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Abstract 

Objective: Cyclooxygenase-2 (COX-2) is an inducible enzyme converting arachidonic acid to prostaglandins and playing 
important roles in inflammatory diseases as well as tumor development. Previous studies investigating the association 
between COX-2 polymorphisms and colorectal cancer (CRC) risk reported conflicting results. We performed a meta-analysis 
of all available studies to explore this association. 

Methods: All studies published up to October 2013 on the association between COX-2 polymorphisms and CRC risk were 
identified by searching electronic databases PubMed, EMBASE, and Cochrane library. The association between COX-2 
polymorphisms and CRC risk was assessed by odds ratios (ORs) together with their 95% confidence intervals (CIs). 

Resu/ts:Jen studies with 6,774 cases and 9,772 controls were included for —1 195A>G polymorphism, 13 studies including 
6,807 cases and 10,052 controls were available for — 765G>C polymorphism, and 8 studies containing 5,121 cases and 7,487 
controls were included for 8473T>C polymorphism. With respect to — 765G>C polymorphism, we did not find a significant 
association with CRC risk when all eligible studies were pooled into the meta-analysis. However, in subgroup analyses by 
ethnicity and cancer location, with a Bonferroni corrected alpha of 0.05/2, statistical significant increased CRC risk was found 
in the Asian populations (dominant model CC+CG vs. GG: OR =1.399, 95%CI: 1.113-1.760, P = 0.004) and rectum cancer 
patients (CC vs. GG: OR = 2.270, 95%CI: 1.295-3.980, P = 0.004; Recessive model CC vs. CG+GG: OR = 2.269, 95%CI: 1.297- 
3.970, P = 0.004). In subgroup analysis according to source of control, no significant association was detected. With respect 
to — 1195A>G and 8473T>C polymorphisms, no significant association with CRC risk was demonstrated in the overall and 
subgroup analyses. 

Conclusions: The present meta-analysis suggests that the COX-2 — 765G>C polymorphism may be a risk factor for CRC in 
Asians and rectum cancer patients. Further large and well-designed studies are needed to confirm this association. 
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Introduction 

Colorectal cancer (CRC) is the second most commonly 
diagnosed cancer with over 1.2 million new cases and 608,700 
deaths in 2008 [1,2]. The highest incidence rate of CRC is found 
in Australia, Europe, and North America [2]. In addition, the 
incidence rate of CRC is rapidly increasing in a number of 
countries within Eastern Asia, such as China [2]. Though the 
exact mechanism of CRC is still unknown, it has been well 
accepted that smoking, obesity, red meat consumption, and 
excessive alcohol consumption are risk factors for CRC [3,4]. 
However, most individuals exposing to these known risk factors 



never develop CRC while many CRC cases develop among 
individuals without those known risk factors, suggesting that other 
factors such as genetic factors also play an important role in the 
development of CRC. 

Cyclooxygenase-2 (COX-2) is an inducible enzyme that 
converts arachidonic acid to prostaglandins, which are potent 
mediators of inflammation. Through the production of prosta- 
glandins, COX-2 is widely regarded as pro-inflammatory factor 
which can be activated by cytokines, mitogens, and growth factors 
at both the transcriptional and post-transcriptional levels [5]. 
Besides, accumulating evidence shows that COX-2 may play a key 
role in tumorigenesis of a variety of human malignancies by 
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stimulating cell proliferation, inhibiting apoptosis, stimulating 
angiogenesis, and mediating immune suppression [6,7,8,9]. The 
human COX-2 gene, mapped to chromosome Iq25.2-q25.3, is 
7.5 kb in length and contains 10 exons [10]. Several potentially 
functional single-nucleotide polymorphisms (SNP), — 765G>C 
(reference SNP ID, rs20417), -11950A (rs689466), and 
8473T>C (rs5275) in the COX-2 gene have been identified. It 
was reported that these three SNPs modulated the inflammatory 
response through affecting gene transcription and/or mRNA 
stability, and consequently contributed to individual variation in 
susceptibility to cancers [11]. CRC is a typical inflammation- 
related malignancy, the pathological progress of CRC is a chronic 
inflammatory process [12]. Hence, it is biologically reasonable to 
hypothesize a potential relationship between the COX-2 gene 
polymorphisms and CRC risk. 

Over the last two decades, a number of molecular epidemio- 
logical studies have been conducted to investigate the association 
between COX-2 -765G>C, -11950A, and 8473T>C 
polymorphisms and CRC risk, but the results remain controversial 
and inconclusive. With respect to — 765G>C polymorphism, a 
meta-analysis by Cao et al [13]. found that individuals carrying the 
GC+CC genotypes were associated with increased risk of CRC 
among Asians (OR = 1.40, 95%CI: 1.11-1.76), however, they 
failed to include the largest sample study by Markar et al. [14] and 
other eligible studies [15,16], which might make their conclusions 
questionable. With respect to -1195G>A, and 8473T>C 
polymorphisms, to the best of our knowledge, no meta-analyses 
on this issue have ever appeared. To derive a more precise 
estimation of the relationship between COX-2 polymorphisms 
and CRC risk, we conducted a meta-analysis of all available case- 
control studies relating the — 765G>C, — 11950A, and 
8473T>C polymorphisms of the COX-2 gene to the risk of 
developing CRC. 

Materials and Methods 

Search strategy 

We conducted a comprehensive literature search in PubMed, 
Embase, and Cochrane library databases up to October 01, 2013 
using the following search strategy: ("colorectal cancer", "CRC", 
"colon cancer" or "rectum cancer") and ("cyclooxygenase-2", 
"COX-2", or PTGS2) and ("polymorphism", "variation", "mu- 
tation", "genotype", or "genetic polymorphism"). There was no 
restriction on time period, sample size, population, language, or 
type of report. All eligible studies were retrieved and their 
references were checked for other relevant studies. The literature 
retrieval was performed in duplication by two independent 
reviewers (Qiliu Peng and Xue Qjn). When multiple publications 
reported on the same or overlapping data, we chose the most 
recent or largest population. When a study reported the results on 
different subpopulations, we treated it as separate studies in the 
meta-analysis. The study was performed according to the proposal 
of Meta-analysis of Observational Studies in Epidemiology group 
(MOOSE) [17]. 

Selection criteria 

The following criteria were used to include published studies: (1) 
Case-control studies which evaluated the association between 
COX-2 polymorphisms and CRC risk; (2) had an odds ratio (OR) 
with 95% confidence interval (CI) or other available data for 
estimating OR (95% CI); and (3) control population did not 
contain malignant tumor patients. Studies were excluded if one of 
the following existed: (1) the design was based on family or sibling 
pairs; (2) the genotype frequency was not reported; (3) there was 



insufficient information for data extraction; and (4) conference 
abstracts, case reports, editorials, review articles, and letters. 

Data extraction 

Two reviewers (Qiliu Peng and Xianjun Lao) independendy 
reviewed and extracted data from all eligible studies. To ensure the 
accuracy of the extracted information, the two investigators 
checked the data extraction results and reached consensus on all of 
the data extracted. If different results were generated, they would 
check the data again and have a discussion to come to an 
agreement. If these two authors could not reach a consensus, 
another author (Xue Qjn) was consulted to resolve the dispute and 
a final decision was made by the majority of the votes. Data 
extracted from eligible studies included the first author, year of 
publication, country of origin, ethnicity, genotyping method, 
matching criteria, source of control, CRC diagnosis criteria, total 
numbers of cases and controls and genotype frequencies of cases 
and controls. Ethnic backgrounds were categorized as Caucasian, 
and Asian. When a study did not state the ethnic descendent or if it 
was not possible to separate participants according to such 
phenotype, the group reported was termed as "mixed ethnicity". 

Statistical analysis 

The strength of the association between COX-2 polymorphisms 
and CRC risk was measured by odds ratios (ORs) with 95% 
confidence intervals (CIs). The significance of the pooled OR was 
determined by Z test and a p value less than 0.05 was considered 
significant. The association of COX-2 polymorphisms with CRC 
risk was assessed using additive models, recessive model, and 
dominant model. 

Heterogeneity in meta-analysis refers to the variation in study 
outcomes between different studies. We used the Q,test and f 
statistics to assess the statistical heterogeneity among studies 
[18,19]. If the result of the Q, test was P h >Q.\ and /<50%, 
indicating the absence of heterogeneity, then a fixed-effects model 
(the Mantel-Haenszel method) was used to estimate the summary 
ORs [20]; otherwise, the random-effects model (the DerSimonian 
and Laird method) was used [21]. To explore the sources of 
heterogeneity among studies, we performed logistic metaregression 
and subgroup analyses. The following study characteristics were 
included as covariates in the metaregression analysis: ethnicity 
(Caucasians versus Asians), source of controls (Hospital-based 
versus Population-based), genotyping methods (PCR-RFLP versus 
not PCR-RFLP), and CRC confirmation (pathologically or 
histologically confirmed versus other diagnosis criteria). Subgroup 
analyses were conducted by ethnicity, cancer location, source of 
control, and HWE in controls. 

Sensitivity analysis was performed by sequential omission of 
individual studies. For each polymorphism, publication bias was 
evaluated using a funnel plot and Egger's regression asymmetry 
test [22]. If publication bias existed, the Duval and Tweedie non- 
parametric "trim and fill" method was used to adjust for it [23]. 
The Bonferroni correction method was used to adjust for multiple 
comparisons. The distribution of the genotypes in the control 
population was tested for HWE using a goodness-of-fit Chi-square 
test. All analyses were performed using Stata software, version 12.0 
(Stata Corp., College Station, TX). M\p values were two-sided. To 
ensure the reliability and the accuracy of the results, two authors 
entered the data into the statistical software programs indepen- 
dently with the same results. 
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Results 

Study characteristics 

Based on the search criteria, 19 studies relevant to the role of 
COX-2 polymorphisms on CRC susceptibility were identified. 
Five of these articles were excluded: two were based on family or 
sibling pairs [24,25], two did not provide allele or genotyping data 
[26,27], and one was a meta-analysis [13]. Manual search of 
references cited in the published studies did not reveal any 
additional articles. As a result, a total of 1 4 relevant studies met the 
inclusion criteria for the meta-analysis 

[14,15,16,28,29,30,31,32,33,34,35,36,37,38] (Figure SI). Among 
them, two of the eligible studies contained data on two different 
ethnic groups, and we treated them independendy [14,31]. 
Therefore, a total of 16 separate comparisons were finally 
included in the meta-analysis. The main characteristics of the 16 
case-control comparisons are summarized in Table 1. Among 
them, 13 studies including 6,807 cases and 10,052 controls were 
available for — 765G>C polymorphism, 10 studies with 6,774 
cases and 9,772 controls for — 1195A>G polymorphism, and 8 
studies containing 5,121 cases and 7,487 controls for 8473T>C 
polymorphism. The sample size of these studies varied consider- 
ably, ranging from 230 to 4,552 individuals. Of all the eligible 
studies, 9 were conducted in Caucasians and 4 were in Asians for 

— 765G>C polymorphism; 8 were conducted in Caucasians and 2 
were in Asians for — 1195A>G polymorphism; all the 8 studies 
were conducted in Caucasians for 8473T>C polymorphism. 
Seven studies were population-based and 9 were hospital-based 
studies. Nine studies in the present meta-analysis did not provide 
definite criteria for the CRC confirmation. Several genotyping 
methods were used, including PCR-RFLP, TaqMan assay, PCR- 
CTTP, and Pyrosequencing . The genotype distributions 
among the controls in two studies were not consistent with 
HWE for -1195A>G [16,32]. 

Meta-analysis results 

Table 2 lists the main results of the meta-analysis of COX-2 - 
1 195A>G polymorphism and CRC risk. There was no evidence 
of significant association between COX-2 -1195A>G polymor- 
phism and CRC risk when all the eligible studies were pooled into 
the meta-analysis (GG vs. AA: OR = 0.902, 95% CI = 0.717- 
1.136, P = 0.380; AG vs. AA: OR = 0.945, 95% CI = 0.835-1.069, 
P = 0.369; GG+AG vs. AA: OR = 0.940, 95%CI = 0.822-1.074, 
P = 0.361, Figure 1; GG vs. AG+AA: OR = 0.891, 
95%CI = 0.786-1.010, P = 0.072). In subgroup analyses by 
ethnicity, cancer location, source of control, and HWE in controls, 
statistical significant association was also not observed in all 
subgroups. 

Table 3 lists the main results of meta-analysis of COX-2 

— 765G>C polymorphism and CRC risk. When all the eligible 
studies were pooled into the meta-analysis, statistical significant 
increased CRC risk was not observed in all genetic models. In 
subgroup analysis according to source of control, significant 
increased CRC risk was also not detected in hospital-based studies 
and population-base studies. However, in subgroup analyses by 
ethnicity and cancer location, after Bonferroni correction for the 
multiple testing (Bonferroni significance threshold P = 0.05 divided 
by the number of ethnicities (n = 2) or cancer types (n = 2): 
P = 0.025), statistical significant increased CRC risk was found in 
Asian populations (dominant model CC+CG vs. GG: OR = 1 .399, 
95%CI: 1.113-1.760, P = 0.004; Figure 2) and rectum cancer 
patients (CC vs. GG: OR= 2.270, 95%CI: 1.295-3.980, 
P = 0.004; Recessive model CC vs. CG+GG: OR = 2.269, 



95%CI: 1.297-3.970, P = 0.004), but not in Caucasian popula- 
tions and colon cancer patients. 

Table 4 lists the main results of meta-analysis of COX-2 
8473T>C polymorphism and CRC risk. There was no evidence 
of significant association between COX-2 8473T>C polymor- 
phism and CRC risk when all eligible studies were pooled into the 
meta-analysis (CC vs. TT: OR = 0.948, 95%CI: 0.843-1.066, 
P = 0.369; TC vs. TT: OR = 1.008, 95%CI: 0.934-1.088, 
P = 0.841; CC+TC vs. TT: OR = 0.995, 95%CI: 0.926-1.070, 
P = 0.899, Figure 3; CC vs. TC+TT: OR = 0.941, 95%CI: 0.842- 
1.051, P = 0.284). In subgroup analyses by ethnicity, cancer 
location, and source of controls, statistical significant association 
was also not observed in all subgroups. 

Heterogeneity analysis 

For the COX-2 — 1195A>G polymorphism, statistical signifi- 
cant heterogeneity among studies was observed when all eligible 
studies were pooled into the meta-analysis (GG vs. AA: P^ = 0.022; 
AG vs. AA: P h = 0.004; GG+AG vs. AA: P h <0.0001). To explore 
the sources of heterogeneity, we first performed subgroup analyses. 
Subgroup analyses by ethnicity, cancer location, and source of 
controls showed that the heterogeneity was still significant in Asian 
populations, hospital-based studies, colon cancer patients and 
rectum cancer patients. Subsequendy, we performed meta- 
regression analysis to further identify the source of heterogeneity. 
Meta-regression analysis indicated that the HWE in controls was 
the major source which contributed to heterogeneity. When we 
excluded two HWE-violating studies [16,32], the heterogeneity 
disappeared (GG vs. AA: P h = 0.833; AG vs. AA: P h = 0.118; 
GG+AG vs. AA: P h = 0.196). However, the significance of the 
summary ORs for COX-2 — 1 195A>G polymorphism in different 
comparison models were not influenced by omitting the two 
studies (Table 2). 

For the COX-2 -765G>C and 8473T>C polymorphisms, 
statistical significant heterogeneity was not detected in the overall 
populations and subgroup analyses. 

Sensitivity analysis 

Sensitivity analysis was performed by sequential omission of 
individual studies. For all these three polymorphisms 
(-1195A>G, -765G>C, and 8473T>C), the significance of 
pooled ORs under all contrast models in both total population and 
subgroup analyses was not influenced excessively by omitting any 
single study (data were not shown). 

For — 1 195A>G polymorphism, sensitivity analysis was further 
performed by omitting those two studies [16,32] in which 
genotype distribution of — 1195A>G polymorphism in the 
controls were significantly deviated from HWE. The significance 
of pooled ORs in both total population and subgroup analyses was 
not influenced by omitting these two studies. 

Publication bias 

Begg's funnel plot and Egger's test were performed to assess the 
publication bias of literatures in all comparison models. The shape 
of the funnel plot did not reveal any evidence of obvious 
asymmetry. Then, the Egger's test was used to provide statistical 
evidence of funnel plot symmetry. The results still did not suggest 
any evidence of publication bias in — 1 195A>G (P = 0.330 for GG 
vs. AA; P= 0.853 for AG vs. AA; P= 0.312 for recessive model 
GG vs. AG+AA; and P= 0.890 for dominant model GG+AG vs. 
AA), -765G>C (P= 0.332 for CC vs. GG; P= 0.815 for CG vs. 
GG; P= 0.389 for recessive model CC vs. CG+GG; and P= 0.703 
for dominant model CC+CG vs. GG), and 8473T>C (P= 0.376 
for CC vs. TT; P= 0.921 for TC vs. TT; P= 0.423 for recessive 
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Table 2. Meta-analysis of COX-2 -1195A>G polymorphism and CRC risk. 



No. of 

Analysis studies Homozygote (GG vs. AA) 



Heterozygote (AG vs. AA) 



Dominant model (GG+AG vs. Recessive model (GG vs. 
AA) AG+AA) 



OR (95% CI) 



OR (95% CI) 



P/Pu 



OR (95% CI) 



P/P h 



OR (95% CI) 



Overall 

Ethnicity 

Caucasian 

Asian 

Cancer 
location 

Colon 

Rectum 

Source of 
control 

HB 

PB 

HWE in 
controls 

Yes 

No 



10 



0.902(0.717-1.136) 0.380/0.022 0.945(0.835-1.069) 0.369/0.004 0.940(0.822-1.074) 0.361 /<0.001 0.891(0.786-1.010) 0.072/0.401 



0.981(0.815-1.181) 0.843/0.833 0.991(0.918-1.069) 
0.809(0.398-1.644) 0.558/0.001 0.881(0.546-1.423) 



0.786(0.537-1.152) 0.217/0.011 0.960(0.751-1.227) 
1.009(0.649-1.568) 0.968/0.010 0.907(0.800-1.029) 



0.918(0.543-1.554) 0.751/0.005 0.975(0.711-1.338) 
0.969(0.799-1.175) 0.750/0.724 0.986(0.910-1.067) 



0.981(0.815-1.181) 0.843/0.833 0.991(0.918-1.069) 
0.809(0.398-1.644) 0.558/0.001 0.881(0.546-1.423) 



0.806/0.118 0.971(0.871-1.083) 0.598/0.096 0.986(0.820-1.185) 0.881/0.909 
0.606/0.005 0.862(0.501-1.484) 0.593/0.001 0.868(0.586-1.286) 0.480/0.034 



0.743/-C0.001 0.928(0.714-1.206) 0.578/-C0.001 0.986(0.655-1.342) 0.339/0.312 
0.130/0.188 0.923(0.763-1.118) 0.413/0.062 1.052(0.747-1.481) 0.771/0.041 



0.877/0.004 0.974(0.685-1.386) 0.449/-C0.001 0.933(0.705-1.384) 0.331/0.132 
0.722/0.149 0.984(0.912-1.062) 0.682/0.130 0.974(0.805-1.179) 0.790/0.804 



0.806/0.118 0.971(0.871-1.083) 0.598/0.196 0.986(0.820-1.185) 0.881/0.909 
0.606/0.005 0.862(0.501-1.484) 0.593/0.001 0.868(0.586-1.286) 0.480/0.034 



P = P values for Z test. P h = P values of Q-test for heterogeneity test. OR, odds ratio; CI, confidence intervals; HB, Hospital-based studies; PB, Population-based studies; 
HWE, Hardy-Weinberg equilibrium. 
doi:1 0.1 371 /journal.pone.0094790.t002 



Study 
ID 



% 

OR (95% CI) Weight 



Caucasian 
Siezenl (2006) 
Siezen2 (2006) 
Hoff (2009) 
Thompson (2009) 
Andersen (2009) 
Pereira (2010) 
Makar(DALS group) (2013) 
Makar(CCFR group) (2013) 
Subtotal (l-squared = 42.4%. p = 0 096) 

Asian 
Tan (2007) 
Li (2013) 

Subtotal (l-squared = 90 8%. p = 0 09^fs 



Overall (l-squared = 70 2%. p = 0.000) 
NOTE: Weights are from random effects analysis 



T 

.422 




0.89 (0.6Z 
0.84 (0.66, 
0.90 (0.66, 
0.87 (0.66, 
0.96 (0.74, 
1.50 (0.95, 
1.13(1.00, 
0.92 (0.80. 
0.97 (0.87. 



1.28)7.40 
1.08)10.15 

1.23) 8.50 
1.14)9.52 

1.24) 9.84 
2.37)5.56 
1.28)13.92 
1.05)13.49 
1.08)78.38 



0.66 (0.55, 0.79)12.13 
1.15(0.87, 1.51)9.50 
0.86 (0.50, 1.48)21.62 

0.94 (0.82, 1.07)100.00 



2 37 



Figure 1. Forest plot of the COX-2 — 1 195G>A polymorphism and CRC risk using a random-effect model (dominant model GG+AG 
vs. AA). 

doi:1 0.1 371 /journal.pone.0094790.g001 
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Table 3. Meta-analysis of COX-2 -765G>C polymorphism and CRC risk. 



Analysis 


No. of 
studies 


Homozygote (CC vs. GG) 


Heterozygote (CG vs. GG) 


Dominant model (CC+CG vs. 
GG) 


Recessive model (CC vs. 
CG+GG) 






OR (95% CI) 




OR (95% CI) 


P/Ph 


OR (95% CI) 


P/P* 


OR (95% CI) 


P/P h 


Overall 


13 


1.197(0.956-1.500) 


0.1 1 7/0.984 


1.048(0.921-1.193) 


0.477/0.031 


1.072(0.958-1.200) 


0.225/0.073 


1.182(0.945-1.478) 


0.143/0.934 


Ethnicity 


Caucasian 


9 


1.196(0.954-1.500) 


0.120/0.968 


1.012(0.934-1.097) 


0.767/0.168 


1.027(0.950-1.111) 


0.504/0.306 


1.180(0.943-1.477) 


0.147/0.890 


Asian 


4 


1.420(0.788-12.931) 


0.805/— 


1.284(0.784-2.104) 


0.320/0.131 


1.399(1.113-1.760) 


O.OO4/0.252 1 .456(0.090-23.477) 


0.791/— 


Cancer 
location 


Colon 


5 


0.983(0.720-1.342) 


0.913/0.889 


1.093(0.982-1.217) 


0.104/0.111 


1.171(0.991-1.382) 


0.063/0.085 


0.966(0.709-1.317) 


0.828/0.856 


Rectum 


4 


2.270(1.295-3.980) 


O.OO4/0.457 1.165(0.788-1.722) 


0.445/0.004 


1.169(0.859-1.591) 


0.321/0.017 
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0.004/0.588 


Source of 
control 


HB 
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P = P values for Z test. P h = P values of Q-test for heterogeneity test. OR, odds ratio; CI, confidence intervals; HB, Hospital-based studies; PB, Population-based studies. 
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model CC vs. TC+TT; and P=0.518 for dominant model Discussion 
CC+TC vs. TT) polymorphisms; Figure 4. 

The present meta-analysis, including 18,702 cases and 27,311 
controls from 16 case-control studies, was conducted systemati- 
cally to evaluate the association between the genetic variants in the 
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Figure 2. Forest plot of the COX-2 — 765G>C polymorphism and CRC risk using a fixed-effect model (dominant model CC+CG vs. 
GG). 
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Table 4. Meta-analysis of COX-2 8473T>C polymorphism and CRC risk. 
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P = P values for Z test. Ph = P values of Q-test for heterogeneity test. OR, odds ratio; CI, confidence intervals; HB, Hospital-based studies; PB, Population-based studies. 
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COX-2 gene and CRC risk. To our knowledge, this is the most 
comprehensive meta-analysis to date to evaluate the association 
between COX-2 polymorphisms and CRC risk. Our results 
showed that the COX-2 — 765G>C polymorphism was associated 
with an increased CRC risk among Asians (dominant model 
CC+CG vs. GG: OR = 1.399, 95%CI: 1.113-1.760, P = 0.004), 
which was in accordance with the previously published meta- 
analysis by Cao et al. [13]. However, with respect to — 1 195A>G 
and 8473T>C polymorphisms, no significant association with 
CRC risk was demonstrated in the overall and subgroup analyses. 



This finding may be biologically plausible. Cyclooxygenases are 
central enzymes in the prostaglandin pathway that convert free 
arachidonic acid into the intermediate prostaglandin H2 which is 
the precursor of prostaglandins, prostacyclin, and thromboxanes. 
Currently, three COX isoenzymes were reported: COX-1, COX- 
2, and COX-3 [39]. COX-2 is normally absent in most cells and 
tissues. It was induced in response to inflammatory cytokines, 
hypoxia, mitogens, hormones, angiogenic growth factors, and 
tumor promoters [40]. COX-2-derived prostaglandins, prostacy- 
clin, and thromboxanes participate in many biologic processes 
such as apoptosis inhibition, inflammation, immune response 



Study 
ID 



% 



OR (95% CI) Weight 



Cox (2004) 
Siezenl (2006) 
Siezen2 (2006) 
Thompson (2009) 
Andersen (2009) 
Pereira (2010) — 
Makar(DALS group) (2013) 
Makar(CCFR group) (2013) 



Overall (l-squared = 0 0%, p = 0.988) 



<> 



0.93 (0.67, 1.30)4.89 
1.02 (0.72, 1.43)4.42 

1.00 (0.79, 1.27)9.11 
1.14 (0.87, 1.48)6.97 

1.01 (0.78, 1.30)7.96 
0.97 (0.62, 1.50)2.72 
0.97 (0.87,1.10)37.46 
0.99 (0.86, 1.14)26.48 
1.00 (0.93,1.07)100.00 



.621 



1.61 



Figure 3. Forest plot of the COX-2 8473T>C polymorphism and CRC risk using a fixed-effect model (dominant model CC+CT vs. TT). 
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Figure 4. Funnel plot analysis to detect publication bias. Each point represents a separate study for the indicated association. A Funnel plot 
for dominant model GG+AG vs. AA of COX-2 -1 195 G>A polymorphism in overall analysis (P = 0.890); B Funnel plot for dominant model CC+CG vs. 
GG of COX-2 -765G>C polymorphism in overall analysis (P = 0.703); C Funnel plot for dominant model CC+CT vs. TT of COX-2 8473T>C 
polymorphism in overall analysis (P = 0.518). 
doi:1 0.1 371 /journal.pone.0094790.g004 



suppression, tumor cell invasion, metastasis, and angiogenesis, 
which are all crucial in the development and progression of cancer 
[39,41,42]. It was shown that polymorphisms in the promoter of 
COX-2 may exert profound effects on gene transcriptional activity 
by altering the binding capacity of certain nuclear proteins, there- 
by affecting expression of COX-2 enzyme [43]. COX-2 

— 765G>C is a functional polymorphism located at 765 bp 
upstream (— 765 bp) from the transcription starting site. It changes 
a putative stimulatory protein (Spl) binding site in the promoter of 
COX-2 between -766 and -761 bp [44], but it creates an E2 
promoter factor (E2F) binding site, leading to high transcription 
activity and increased COX-2 expressions which might be 
involved in the development of cancers [45] . More importantly, 
the homozygous variant genotype COX-2 — 765CC has been 
shown associated with increased risk for many different types of 
cancers, including breast cancer [46], ovarian cancer [47], 
hepatocellular carcinoma [48] and lung cancer [49]. 

In subgroup analysis by ethnicity, the COX-2 — 765G>C 
polymorphism presented a risk factor for CRC in Asian 
populations, but not in Caucasians. The inconsistent data among 
the different ethnicities may indicate different effects of the COX- 
2 — 765G>C polymorphism on CRC risk in different ethnic 
genetic backgrounds. Nevertheless, owing to the limited number of 
relevant studies among Asian population included in this meta- 
analysis, the observed positive association between COX-2 

— 765G>C polymorphism and CRC risk in Asians is likely to 
be caused by chance because study with small sample sizes may 
have insufficient statistical power to detect a slight effect or may 
have generated a fluctuated risk estimate. Currently there are only 
4 studies on COX-2 — 765G>C polymorphism and CRC risk 
among Asian population [28,30,32,33]. Therefore, the positive 
results of the Asain population should be interpreted with caution. 

Studies have suggested that cancers of the rectum and colon 
might be distinct tumors because they have a differing prevalence 
with a difference in clinical presentation, prognosis and possibly in 
genetic and environmental epidemiology. Thus, the COX-2 

— 765G>C polymorphism might influence carcinogenesis of 
colorectal tissues in a site-specific manner. Therefore, we carried 
out subgroup analysis according to cancer location. Our results 

References 

1. Akin O, Brcnnan SB, Dcrshaw DD, Ginsberg MS, Gollub MJ, ct al. (2012) 
Advances in oncologic imaging: update on 5 common cancers. CA Cancer J Clin 
62: 364-393. 



suggested a significant increased CRC risk in rectum cancer 
patients (CC vs. GG: OR= 2.270, 95%CI: 1.295-3.980, 
P = 0.004; Recessive model CC vs. CG+GG: OR = 2.269, 
95%CI: 1.297-3.970, P = 0.004) but not in colon cancer subjects, 
which was consistent with the results of the large sample study by 
Markar et al. [14]. Our findings add further data to evidence that 
colon and rectal cancers have different etiologies. 

In our meta-analysis, several limitations should be acknowl- 
edged. First, in subgroup analyses by ethnicity and cancer location, 
the sample size of population was relatively small for subgroup 
analyses, which may lead to relatively weak power to detect the 
real relationship. Second, our results were based on unadjusted 
estimates. We did not perform the analysis adjusted for other 
covariates such as smoking, drinking, obesity, red meat consump- 
tion, and so on, because of the unavailable original data of the 
eligible studies. 

In conclusion, our meta-analysis provided a more precise 
estimation based on larger sample size compared with the 
individual studies and previous meta-analysis. Our study suggested 
that COX-2 — 765G>C polymorphism might contribute to 
colorectal cancer risk, especially in Asian populations and the 
rectum cancer patients. In order to further verify our findings, 
large well designed epidemiological studies are warranted. 

Supporting Information 

Figure SI Flow diagram of included studies for this 
meta-analysis . 

(TIF) 

Checklist SI 

(DOC) 

Author Contributions 

Conceived and designed the experiments: QP XL WT SL XQ. Performed 
the experiments: HL JW SY JS. Analyzed the data: ZC XQ QP SY. 
Contributed reagents/materials/analysis tools: QP XQ. Wrote the paper: 
QP XL SY. 



2. Jcmal A, Bray F, Center MM, Ferlay J, Ward E, et al. (2011) Global cancer 
statistics. CA Cancer J Clin 61: 69-90. 



PLOS ONE | www.plosone.org 



8 



April 2014 | Volume 9 | Issue 4 | e94790 



COX-2 Polymorphisms and Colorectal Cancer 



3. Theodoratou E, Montazeri Z, Hawkcn S, Allum GC, Gong J, ct al. (2012) 
Systematic meta-analyses and field synopsis of genetic association studies in 
colorectal cancer. J Natl Cancer Inst 104: 1433-1457. 

4. Burn J, Mathers J, Bishop DT (2013) Genetics, inheritance and strategics for 
prevention in populations at high risk of colorectal cancer (CRC). Recent Results 
Cancer Res 191: 157-183. 

5. Hla T, Bishop-Bailey D, Liu CH, Schaefers HJ, Trifan OC (1999) 
Cyclooxygenase-1 and -2 isoenzymes. IntJ Biochcm Cell Biol 31: 551-557. 

6. Kirschcnbaum A, Liu X, Yao S, Levine AC (2001) The role of cyclooxygenase-2 
in prostate cancer. Urology 58: 127—131. 

7. Fujita H, Koshida K, Keller ET, Takahashi Y, Yoshimito T, et al. (2002) 
Cyclooxygenasc-2 promotes prostate cancer progression. Prostate 53: 232-240. 

8. Nithipatikom K, Isbell MA, Lindholm PF, Kajdacsy-Balla A, Kaul S, et al. 
(2002) Requirement of cyclooxygenase-2 expression and prostaglandins for 
human prostate cancer cell invasion. Clin Exp Metastasis 19: 593—601. 

9. Wang W, Bergh A, Dambcr JE (2005) Cyclooxygenase-2 expression correlates 
with local chronic inflammation and tumor neovascularization in human 
prostate cancer. Clin Cancer Res 11: 3250—3256. 

10. Tazawa R, Xu XM, Wu KK, Wang LH (1994) Characterization of the genomic 
structure, chromosomal location and promoter of human prostaglandin H 
synthase-2 gene. Biochcm Biophys Res Commun 203: 190-199. 

11. Dong J, Dai J, Zhang M, Hu Z, Shen H (2010) Potentially functional COX-2- 
1195G>A polymorphism increases the risk of digestive system cancers: a meta- 
analysis. J Gastroenterol Hepatol 25: 1042—1050. 

12. Itzkowitz SH, Yio X (2004) Inflammation and cancer IV. Colorectal cancer in 
inflammatory bowel disease: the role of inflammation. AmJ Physiol Gastrointcst 
Liver Physiol 287: G7-17. 

13. Cao H, Xu Z, Long H, Li XQ, Li SL (2010) The -765C allele of the 
cyclooxygcnase-2 gene as a potential risk factor of colorectal cancer: a meta- 
analysis. Tohoku J Exp Med 222: 15-21. 

14. Makar KW, Poole EM, Reslcr AJ, Seufert B, Curtin K, et al. (2013) COX-1 
(PTGS1) and COX-2 (PTGS2) polymorphisms, NSAID interactions, and risk of 
colon and rectal cancers in two independent populations. Cancer Causes 
Control. 

15. Daraei A, Salchi R, Mohamadhashem F (2012) PTGS2 (COX2) -765G>C gene 
polymorphism and risk of sporadic colorectal cancer in Iranian population. Mol 
Biol Rep 39: 5219-5224. 

16. Li S, Zhao X, Wu Z, Li Y, Zhu L, ct al. (2013) Polymorphisms in arachidonic 
acid metabolism-related genes and the risk and prognosis of colorectal cancer. 
Fam Cancer. 

17. Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, et al. (2000) Meta- 
analysis of observational studies in epidemiology: a proposal for reporting. Meta- 
analysis Of Observational Studies in Epidemiology (MOOSE) group. JAMA 
283: 2008-2012. 

18. Higgins JP, Thompson SG (2002) Quantifying heterogeneity in a meta-analysis. 
Stat Med 21: 1539-1558. 

19. Higgins JP, Thompson SG, Decks JJ, Altman DG (2003) Measuring 
inconsistency in meta-analyses. BMJ 327: 557—560. 

20. Mantel N, Haenszcl W (1959) Statistical aspects of the analysis of data from 
retrospective studies of disease. J Natl Cancer Inst 22: 719—748. 

21. DcrSimonian R, Laird N (1986) Meta-analysis in clinical trials. Control Clin 
Trials 7: 177-188. 

22. Egger M, Davey Smith G, Schneider M, Minder C (1997) Bias in meta-analysis 
detected by a simple, graphical test. BMJ 315: 629-634. 

23. Duval S, Tweedie R (2000) Trim and fill: A simple funnel-plot-bascd method of 
testing and adjusting for publication bias in meta-analysis. Biometrics 56: 455- 
463. 

24. Ross J, Lockctt L, Brookes D, Tabor B, Duesing K, et al. (2013) An association 
between the PTGS2 rs5275 polymorphism and colorectal cancer risk in families 
with inherited non-syndromie predisposition. Eur J Hum Genet. 

25. Wang J, Joshi AD, Corral R, Siegmund KD, Marchand LL, ct al. (2012) 
Carcinogen metabolism genes, red meat and poultry intake, and colorectal 
cancer risk. IntJ Cancer 130: 1898-1907. 

26. Kwagyan J, Apprey V, Ashktorab H (20 1 2) Linkage disequilibrium and 
haplotypc analysis of COX-2 and risk of colorectal adenoma development. Clin 
Transl Sci 5: 60-64. 

27. Goodman JE, Bowman ED, Chanock SJ, Alberg AJ, Harris CC (2004) 
Arachidonate lipoxygenase (ALOX) and cyclooxygenase (COX) polymorphisms 
and colon cancer risk. Carcinogenesis 25: 2467—2472. 

28. Hamajima N, Takezaki T, Matsuo K, Saito T, Inoue M, et al. (2001) Genotype 
Frequencies of Cyclooxygeneasc 2 (COX2) Rare Polymorphisms for Japanese 
with and without Colorectal Cancer. Asian PacJ Cancer Prev 2: 57—62. 



29. Cox DG, Pontes C, Guino E, Navarro M, Osorio A, et al. (2004) Polymorphisms 
in prostaglandin synthase 2/cyelooxygenasc 2 (PTGS2/COX2) and risk of 
colorectal cancer. Br J Cancer 91: 339—343. 

30. Koh WP, Yuan JM, van den Berg D, Lec HP, Yu MC (2004) Interaction 
between cyclooxygenase-2 gene polymorphism and dietary n-6 polyunsaturated 
fatty acids on colon cancer risk: the Singapore Chinese Health Study. BrJ Cancer 
90: 1760-1764. 

31. Siezen CL, Bucno-dc-Mcsquita HB, Peeters PH, Kram NR, van Doeselaar M, 
ct al. (2006) Polymorphisms in the genes involved in the arachidonic acid- 
pathway, fish consumption and the risk of colorectal cancer. IntJ Cancer 119: 
297-303. 

32. Tan W, Wu J, Zhang X, Guo Y, Liu J, ct al. (2007) Associations of functional 
polymorphisms in cyclooxygenase-2 and platelet 12-lipoxygenase with risk of 
occurrence and advanced disease status of colorectal cancer. Carcinogenesis 28: 
1197-1201. 

33. Xing LL, Wang ZN, Jiang L, Zhang Y, Xu YY, et al. (2008) Cyclooxygenase 2 
polymorphism and colorectal cancer: -765G>C variant modifies risk associated 
with smoking and body mass index. WorldJ Gastroenterol 14: 1785—1789. 

34. Hoff JH, tc Morsche RH, Roelofs HM, van der Logt EM, Nagengast FM, et al. 
(2009) COX-2 polymorphisms -765G— >C and -1195A— >G and colorectal 
cancer risk. WorldJ Gastroenterol 15: 4561—4565. 

35. Thompson CL, Plummer SJ, Merkulova A, Cheng I, Tucker TC, et al. (2009) 
No association between cyclooxygenasc-2 and uridine diphosphate glucurono- 
syltransfcrase 1A6 genetic polymorphisms and colon cancer risk. 
WorldJ Gastroenterol 15: 2240-2244. 

36. Andersen V, Ostergaard M, Christensen J, Ovcrvad K, Tjonneland A, et al. 
(2009) Polymorphisms in the xenobiotic transporter Multidrug Resistance 1 
(MDR1) and interaction with meat intake in relation to risk of colorectal cancer 
in a Danish prospective case-cohort study. BMC Cancer 9: 407. 

37. Iglesias D, Ncjda N, Azcoita MM, Schwartz S, Jr., Gonzalcz-Aguilcra JJ, et al. 

(2009) Effect of COX2 -765G>C and c.3618A>G polymorphisms on the risk 
and survival of sporadic colorectal cancer. Cancer Causes Control 20: 1421- 
1429. 

38. Percira C, Pimentel-Nunes P, Brandao C, Moreira-Dias L, Mcdeiros R, ct al. 

(2010) COX-2 polymorphisms and colorectal cancer risk: a strategy for 
chemoprevention. Eur J Gastroenterol Hepatol 22: 607-613. 

39. Mullcr-Dcckcr K (2011) Cyclooxygenase-dependent signaling is causally linked 
to non-melanoma skin carcinogenesis: pharmacological, genetic, and clinical 
evidence. Cancer Metastasis Rev 30: 343-361. 

40. Miyashita M, Ito T, Sakaki M, Kajiwara A, Nozawa H, et al. (2012) Genetic 
polymorphism in cyclooxygenase-2 promoter affects hepatic inflammation and 
fibrosis in patients with chronic hepatitis C.J Viral Hepat 19: 608—614. 

41. Rouzcr CA, Marnett LJ (2011) Endoeannabinoid oxygenation by eyclooxygen- 
ases, lipoxygenases, and cytochromes P450: cross-talk between the eicosanoid 
and endoeannabinoid signaling pathways. Chcm Rev 111: 5899-592 1 . 

42. Wang YX, GaoJX, Wang XY, Zhang L, Liu CM (2012) Antiproliferative effects 
of selective cyclooxygenase- 2 inhibitor modulated by nimotuzumab in estrogen- 
dependent breast cancer cells. Tumour Biol 33: 957—966. 

43. Brosens LA, Iacobuzio-Donahue CA, Keller JJ, Hustinx SR, Carvalho R, et al. 
(2005) Increased cyclooxygcnase-2 expression in duodenal compared with 
colonic tissues in familial adenomatous polyposis and relationship to the -765G - 
> C COX-2 polymorphism. Clin Cancer Res 11: 4090^096. 

44. Papafili A, Hill MR, Brull DJ, McAnulty RJ, Marshall RP, ct al. (2002) 
Common promoter variant in cyclooxygcnase-2 represses gene expression: 
evidence of role in acute-phase inflammatory response. Arteriosclcr Thromb 
Vase Biol 22: 1631-1636. 

45. Szczcklik W, Sanak M, Szczcklik A (2004) Functional effects and gender 
association of COX-2 gene polymorphism G-765C in bronchial asthma. 
J Allergy Clin Immunol 114: 248-253. 

46. Gao J, Ke Q, Ma HX, Wang Y, Zhou Y, et al. (2007) Functional polymorphisms 
in the cyclooxygenase 2 (COX-2) gene and risk of breast cancer in a Chinese 
population. J Toxicol Environ Health A 70: 908-915. 

47. Agachan Cakmakoglu B, Attar R, Kahraman OT, Dalan AB, Iyibozkurt AC, et 
al. (2011) Cyclooxygenase-2 gene and epithelial ovarian carcinoma risk. Mol 
Biol Rep 38: 3481-3486. 

48. Akkiz H, Bayram S, Bekar A, Akgollu E, Ulgcr Y (2011) Functional 
polymorphisms of cyclooxygenase -2 gene and risk for hepatocellular carcinoma. 
Mol Cell Biochcm 347: 201-208. 

49. Coskunpinar E, Eraltan IY, Turna A, Agachan B (201 1) Cyclooxygenase- 2 gene 
and lung carcinoma risk. Med Oncol 28: 1436-1440. 



PLOS ONE | www.plosone.org 



9 



April 2014 | Volume 9 | Issue 4 | e94790 



